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ABSTRACT 

We have taken a deep, moderate-resolution Keck/Deimos spectra of QSO, CFHQS2329, at 
z=6.4. At the wavelength of Lya, the spectrum shows a spatially-extended component, which 
is significantly more extended than a stellar spectrum, and also a continuum part of the spec- 
trum. The restframe line width of the extended component is 21 ±7 (A), and thus smaller than 
that of QSO (52±4A) , where they should be identical if the light is incomplete subtraction 
of the QSO component. Therefore, these comparisons argue for the detection of a spatially 
extended Lya nebulae around this QSO. This is the first z>6 QSO that an extended Lya halo 
has been observed around. Careful subtraction of the central QSO spectrum reveals a lower 
limit to the Lya luminosity of (1.7±0.1)xl0'*'^ erg s^^. This emission may be from the the- 
oretically predicted infalling gas in the process of forming a primordial galaxy that is ionized 
by a central QSO. On the other hand, if it is photoionized by the host galaxy, an estimated 
star-formation rate of >3.0 Mq yr~^ is required. 

If we assume the gas is virialized, we obtain dynamical mass estimate of Mdyn=^-'2 
xlO^^M0. The derived MBnlMhost is 2.1 xlO~^, which is two orders smaller than those 
from more massive z^6 QSOs, and places this galaxy in accordance with the local M-cr re- 
lation, in contrast to a previous claim on the evolution of M-tr relation at z~6. We do not 
claim evolution or non-evolution of the M-cr relation based on a single object, but our result 
highlights the importance of investigating fainter QSOs at z^6. 
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1 INTRODUCTION 

Understanding the first stage of the formation of galaxies is one 
of the central issues in observational astronomy. Galaxy formation 
models predict that an early stage inevitably involves a spatially ex- 
tended distribution of infalling, cold gas (Haiman & Rees 2001). If 
such gas is ionized by a central luminous quasar (QSO), or an initial 
starburst, such gas should be observed as extended Lya emission 
in the high-z Universe. 

However, due to instrumental limitations, extensive searches 
for extended Lya emissi on or Lyg blobs were mostly conducted 



at at zr^2-3 in the past jSteidel et al.ll200d iMatsudaet al.|l2004 
iNilsson et aUbOOd Ismith & Ja rvis"2007: Smith et al. "2009'), with 
only a few exceptio ns at z=6.595 (Ouchi et al. 2009) and at z=4.5 
( lBunkeretal]|2003l) . It is worth noting that some of these de- 
tections are associated with a bright Q SO dBunker et alj l2003l : 
IWeidingeret"Zll2005l : iBarrio et aljboOSh . If a QSO is a heating 
source, it is also a subject of interest to understand what role a cen- 



tral black hole plays in such an early stage of the galaxy formation, 
perhaps leading to the ti ght correlations obser ved at low redshift 
such as the M-cr relation (Magorrian et al. "1998'), and the starburst- 
AGN connection tMaiolino et al. 1997 : Goto 2005. 2006) . 



However, at z=3, the age of the Universe was already 2.2 Gyrs, 
and thus, it may be too late to search for a primordial stage of a 
galaxy formation (Scannapieco et al. 2003). In Goto et al. (200^, 
we found a spatially-extended structure around a QSO at z=6.4. 
This may be the first example of Lya blob around a luminous 
QSO at z>6. However, the detection was in the z'-band image 
of Subaru/Suprime-Cam. Therefore, it was not clear if the ex- 
tended structure was Lya emission, or continuum emission from 
the QSO host galaxy. In this work, we performed a deep, moderate- 
resolution spectroscopy with the Keck/Deimos to clearly separate 
the two cases. Unless otherwise stated, we adopt the WMAP cos- 
mology: {h,nm,^L) = (0.7,0.3,0.7). 



2 Goto 

Table 1. Target information adopted from lWillott et alj j2007l) and lGoto et all j2009 l). 
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J 


CFHQS J232908.28-030158.8 


6.417±0.002 


25.54±0.02 21.165±0.003 21.683±0.007 


21.56±0.25 
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Figure 1. Subtracting stellar spectrum from QSO spectrum. The top panel 
shows the QSO spectrum. The middle panel is a reference PSF stellar spec- 
trum. The bottom panel shows the residuals from the subtraction of the 
PSF spectrum and the smooth extended component from the QSO spec- 
tram. Pixel scale in spatial direction is 0.1185" pix^^. The red bar shows a 
scale of 2". 



2 OBSERVATION 



Our target i s QSO CFHQS J23 29-0301(Tablem lw"ilIott et al.l200"7h 
atz=6.417 jWillott et al.l2010h . This QSO is kn own to be in a dens e 
environment surrounded by 7 LBG candidatesC Utsumi et al.l201oh . 

We obtained moderate resolution spectra of the QSO using 
Keck/Deimos on the nights of September 12 and 13, in 2010. 
The details of the observation and data reduction are described in 
lGotoetalj|20Tlh . Briefly, we used the 830 lines mm ^ grating and 
the OG550 order cut filter with the central wavelength of 8500A . 
The slit width was 1.0" with 0.47 A pixel" ^, giving a resolving 
power of R ~3600. The position angle of the slit was -30 deg. The 
wavelength coverage was 6000 A to lOOOOA. The spatial resolution 
is 0.1185" pixel"^. The total integration time was 5.5 hours. 

We used the Deep2 pipeline to reduce the data, except the 
background subtraction, which we did manually to remove ghost 
features of the 830 grating carefully. Wavelength calibration is 
based on the HeNeAr lamp. 

For flux calibration, the spectrophotometric standard G191- 
B2B was observed and used to correct the spectral shape. Absolute 
flux calibration was achieved by passing the spectra of the QSO 
through the Subaru z' filter and normalizing to match the observed 
z' magnitude in Table[T] We independently performed the absolute 
flux calibration using the standard G191-B2B, obtaining only 4% 
larger absolute flux. 



3 ANALYSIS 

3. 1 Removal of the QSO PSF 

To investigate the extended component of the spectrum, we need to 
subtract the central point source spectrum, which is often brighter 
by a factor of >10 for a QSO. We chose a spectrum of a bright 
star in the same Deimos mask as the QSO and which was rela- 
tively free of the ghost features of the 830G grism; this spectrum is 
shown in the middle panel of Fig[T] The PSF of the QSO and the 
extended background were decomposed (in the spatial direction) 



Table 2. Properties of Lya emission from the extended region. 



FWHM (A) Flux (erg cm ^) Luminosity (erg s ^) Velocity Disp. (km s ^) 



21±7 (3.6±0.2) xlO-i'^ (1.7±0.1)xl043 



301±99 



using the IRAF code specinholucy based on a two-channel restora- 
tion technique ( Lucy & W alsh 2003 ). The technique restores the 
spatial profile of the PSF and a smooth background wavelength- 
by-wavelength. It is based on a two-channel restoration algorithm 
that restores a point spread function (PSF)-like component in a 2D 
spectrum and an underlying extended background component. It 
has already been successfully used to subtract point-source spec- 
tra in highly inhomogeneous b ackgrounds, such as high-z SNe la 
embedded in their host galaxy jBlondin et alj|200^ . Therefore, it 
is also suitable for our case of separating QSO from the extended 
structure. 

The top panel of Fig[T] shows the reduced 2D spectrum of 
J2329. The middle panel shows the spectrum of the reference PSF 
star. The PSF size (FWHM) measured from a star in a slit is 0.65". 

The restored background component was smoothed by a 
Gaussian (see ILucv & Walslj 120031. for details) of FWHM 16.5 
pixels (1.95"). On subtracting this smooth component from the 
spectrum (with the QSO PSF removed), structure at spatial scales 
intermediate between that of the PSF and the smooth extended 
structure was revealed. The bottom panel of Fig[T] shows this re- 
vealed component. 

In the bottom panel of Fig[T] at the peak of the QSO spec- 
trum, there remain positive and negative residuals resulting from 
incomplete PSF removal. This has often been seen in PSF removal 
and usually arises when the spatial profile of the PSF star does not 
quite match that of the target. The residuals are at the level of a few 
percent of the peak of the QSO PSF. 

Noteworthy in Fig[T]is the extended flux over the wavelength 
range 9006-9035A that remains from the decomposition. The peak 
of this flux corresponds to Lya wavelength at the redshift of the 
QSO. The red bar of length 2" in the figure shows the spatial 
scale for reference, indicating that the remaining flux is spatially 
extended over an extent as large as ~4". In Fig|2l we show the 
residual extended flux from the two-channel restoration summed in 
the wavelength direction over the extent 9006 to 9035A. The er- 
ror bars include Poisson noise from the QSO and star spectra, and 
the background noise. Compared with the errorbars, remaining flux 
from -1.7 to -1-1.6 arcsec appears to be significant. 

To further verify this finding, in FigjS] we compare the spatial 
flux profile in the wavelength range between 9006 and 9035A of 
the QSO spectrum (corresponding to Lya emission) in the black 
solid line, and that between 9058 and 9087A(i.e., remote from the 
Lyof line) in the red dotted line. Both of these wavelength ranges are 
chosen so that they do not include strong sky emission lines. The 
flux profile of the star in the range of 9006-9035A is also plotted 
in the blue dotted line. Compared with the associated errorbars, 
the flux profile of the QSO at the Lya wavelength is significantly 
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Figure 2. The residuals on the two channel restoration of the QSO 2D 
spectrum are shown summed in wavelength from 9006 to 9035A. The error 
bars include those of the QSO, the star, and the background noise. 
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Figure 3. Comparison of the spatial profiles of the QSO spectrum over 
the Ly-a line between 9006 and 9035A (black solid line), and over the 
continuum between 9058 and 9087A (red dotted line). The blue dotted line 
show the spatial profile of a star between 9006 and 9035A. 



more extended than that of the continuum, and the star. We also 
performed a 2-dimensional subtraction using the QSO continuum 
as a spectral PSF, interpolating the PSF shape from the continuum 
on both sides of the Lya line, obtaining a very similar result as 
in Fig[T] Note that the continuum region may contain an extended 
component from the underlying galaxy. This test was performed in 
case the optical distortions in the spectrograph differ between the 
position of the star used as a PSF and QSO slit position. However, 
the spatial profiles of the star and QSO continuum are very similar 
(see Fig.O, suggesting such effects are small. 

We note that the extended emission may be stronger above 
than below the QSO spectrum in Fig[T] Corresponding asymmetry 
can be found in Fig|2] This can also be seen in Fig|3] where the ex- 
cess from the continuum subtraction is larger at the positive offset. 
Also the emission may be tilted in wavelength-position space as re- 
vealed in the bottom panel of Fig[T] This may be a sign of rotation 
in the extended Lya emitting region, but needs confirmation from 
a better spatial resolution spectrum. 



Figure 4. One-dimensional spectra of the extended region of CFHQSJ2329 
after subtracting the central QSO. The gray line shows spectrum in the orig- 
inal resolution. The black line shows binned spectrum with 10 pixel box. 
The magenta dash-dotted lines show the sky spectrum in an arbitrary unit. 
The blue dotted line shows the arbitrary-scaled spectrum of the QSO. 

3.2 One dimensional spectrum 

On the basis of this evidence, we claim that at least the Lya region 
of the spectrum is more extended than the PSF from the star and 
the continuum region of the QSO spectrum. Next, we extracted a 
one-dimensional spectrum from the 2D image by masking out the 
central 10 pixels (1.1"), where the residuals on the PSF match are 
strong. Note that due to this masking, the computed flux and lumi- 
nosity are lower limits. 

In Fig|4l we show the one-dimensional spectrum. The gray 
line shows the spectrum at the original resolution. The black line 
shows the same spectrum binned with a 10 pixel box. The ma- 
genta line shows the scaled sky spectrum. As expected from Fig[T] 
a strong peak can be seen at around Lya emission. Due to neutral 
hydrogen absorption on the bluer side, the line has clear asymmet- 
ric profile with an extended red wing that is typical for a high-z 
Lya line. In the blue dotted line, we show the arbitrary-scaled QSO 
spectrum before the PSF subtraction. The Lya line profile of the 
QSO is much broader. This difference confirms that the residual 
in the bottom panel of Fig[T]is not part of the (unsubtracted) QSO 
spectrum, but indeed spectrum of a different nature. 

The total flux measured between 9006 and 9035A region (ex- 
cept central 1") is (3.6±0.2) xl0~^^ erg cm~^s~^ as shown in 
Table |2] Compared with errors measured from 9058-9087A region 
of the spectrum, this is a 17 a detection. There is no significant 
flux detected other than the Lya wavelength. The flux corresponds 
to the lower limit of Lya luminosity of (I.7±0.1)xlQ''^ erg s^^ 
(c.f., the Lya luminosity of the QSO is 6.2x10'*'' erg s"^). This 
is sm aller than photom etrically estimated value of 2.0X10'*'' erg 
s~'n( iGotoet"al]l2009i ). but consistent with each other, since our 
estimate is a lower limit due to the central mask. For comparison, 

^ iGoto et ai] )2009l) had a computational error We quoted a corrected value 
here. 
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Figure 5. Mgu — a relation. The data points except for J2329 are 
adopted from Wang et al. (2010). The diamonds represent the local galax- 
ies ( Tremaine et a l. 2002), whose local Mbh ~ relationship in the black 
dashed line is log (Mbh/Mq) = 8.13 + 4.021og(o-/200 km s^i). The 
blue squares and triangles are for the z~6 and 1.4 z ^ 5 QSO sam- 
ples, respectively, with a derived from th e observed CO line width using 
o- FWHMI2.Z^ JWang et aljr20ldh . The red circle show J2329 at 
z=6.4 measured in this work based on the dispersion of t he Lyg line. The 
purpl e dotted line shows a simulated M-cr relation at z=6 iRobertson et alj 
l200el) . 



-6.5 is L* ~3.9±0.2x 10''=' erg 



touchi et al."2009^ and Lya blobs at z=3.1 is 
(Matsuda et al. 2004, 2006). 



1 X erg 



the l uminosity of a Lyg blob at z 

s 
s 

The FWHM estimated through a Gaussian fit is 21 ±7 (A) in 
the restframe. This narrow width of the line also suggests that the 
Lyof emission did not originate from the broad line regions of the 
QSO (FWHM of 52±4A). We do not recognize the presence of a 
clear rotation curve, or multiple velocity components. In terms of 
velocity dispersion in the restframe, this corresponds to 301 ±99 
km s~^. For comparison, the median velocity d ispersion of Lya 
blobs at z=6 .6 is 110 km s~^ Jouchi et al.1120091) , and at z=3.1 is 



330 km s^^ ('Matsu da et alj2004 l200^ 



If the nebula forms a single virialized system with a veloc- 
ity dispersion of 301 km s~^ in 2" radius (11 kpc at z=6.4), we 
estimate a virial mass of 1 .2 x lO^^M© (5/3x3cr^7?/G). For com- 
parison, iGoto et al.l l2009h estimated a stellar mass from 6.2 x 10 
to 1 . 1 X 10^" Af0 depending on the star-formation history. If the ve- 
locity dispersion of the nebula reflects the dynamics of the host 
galaxy, it suggests that a massive galaxy is already harboring a QSO 
at z=6.4. 



3.3 Non-detection of Continuum 

iGoto et alj j2009l) detected 2.5 x 10"^'-* erg cm"^ s"^ A"^ of flux 
density from the extended component in Zr filter (Ae=988nm) at 
3 (J. We did not detect any continuum from the extended compo- 
nent other than the Lya emission. However, the noise level of our 
continuum is 5.5xlO~^' erg cm~^ s"'^ A~^ pixel"^. Therefore, 
we cannot indicate the detection or non-detection of the possible 
continuum emission from the extended component. 

Since we did not detect continuum from the extended compo- 
nent, we cannot obtain an equivalent width (EW) solely from the 
spectra. However, if we use the continuum level of Zr band, we ob- 



tain the restframe EW of 19±6 A. This is not a partic ularly large 
value for Lya blobs. For example, all 18 Lya blobs in lSaito et al.l 
|2008) had EWrest > lOOA (but not all of them are associated 
with a QSO). 



4 DISCUSSION 

4.1 Comparison to Lya blobs around low-z QSOs 

It is informative to compare the size of the extended Ly-a emis- 
sion around CFHQS J232908.28-030158.8 with th at of similar 
structu res, Ly-a blobs around lower redshift QSOs. iBarrio et al.l 
( l2008l) 's Lya nebula around a z=2.48 QSO had FWHM^est of 
370±63 km/s, with 40 kpc radius of extension (in 3 separate parts 
though). The Lya h alo around a radio-qui e t QSO at z=3.04 also ex- 
tended to ~30 kpc jWeidinger et alj2005l) . rBunker et alj j2003h re- 
ported a 5"-extended, FWHMrest of 181 km/s Lya nebula around 
a QSO at z=4.5. Compared to these, our z=6.4 Lya nebula has a 
much larger FWHMrest of 707±232 km/s, perhaps reflecting a 
larger Mhost- Although the 11 kpc of extension is smaller than 
the measured sizes at lower-z, surface brightness dimming is much 
greater at z=6.4, and thus, the size could be larger than measured. 
There exist only a few examples of extended Lya halos around 
QSOs. It is important to construct a larger sample to correctly un- 
derstand the evolution of the Lya halo around QSOs. 



4.2 What is the physical origin of this Lya nebula? 

Since there exists a QSO at the center of this nebula, an obvi- 
ous ex planation is halo gas photoionized by the QSO. Following 
IYu & L u (2005), we estimate the ionizing photon em ission rate of 
this QSO using absolute magnitude M ° =-25.23 ( ^Willott et al. 

I 1450A ' ' 

120071) . We count photons with energy in the range 13.6-54.4 eV 
as io nizing photons. W e assumed average QSO SEDs at z > 3 
from lTelferet"d] ilOOt). To be conservative, we used the SED of 
radio-loud QSOs, which produce smaller ionizing photon rates than 
radio-quiet QSOs. The rate obtained for the QSO is 1.6x 10'"^''s~\ 
On the other hand, Lya photon rate emitted by the nebula is several 
10^"^ s~^. Therefore, there are orders more ionizing photons pro- 
duced by the QSO to illuminate the Lya nebula. A cold accretion 
of neutral gas tearkan a & Loeb 2 003), or recently fo und re-scatted 
Lya photons by neutral hydrogen faaves et al.l201ll) are other can- 
didates to explain Lya blobs. In our case, however, these scenarios 
are unlikely because neutral hydrogen cannot survive exposure to 
the abundant ionizing photons. 

In addition, star-formation in the host galaxy can also con- 
tribute dTaniguchi & Shiovall2000l : lOhvama et alj|2003l) . The star- 
formation rate (S FR) can be estimated using the following relation 
([Kennicutt 1998 dTaniguchi et alj|2007l) . 



SFR{Lya) = 9.1 x W'^'-^ L{Lya)MQyr- 



(1) 



Based on the Lya luminosity, we estimate the SFR(Lya) is >3.0 
Mq yr"^. This is not a high SFR, and thus, can be expected from 
a young star-forming galaxy. Therefore, it is possible that the star- 
formation in the host galaxy contributes to ionizing the Lya nebula. 

4.3 Mbh — cr relation at z >6 

The virial mass estimate provides an interesting opportunity to in- 
vestigate the MsHlMhost ratio at z> 6. The black hole mass of 
this QSO was measured to be 2.5±0.4 x lO^M© based on the Mgll 
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line and L ° jWillott et alj 20101) . Combined with our virial 

3000A ' ^ 

mass of 1.2 x IO^^Mq, we obtain MsH/Mhost of 2.1 xlO~''.This 
value is much smaller than thos e previously mea sured from brighter 
QSOs at z~6. For example. iWang et alj j2O10l) obtained a median 
MBHiMbuige ratlo of 0.022 for bright {Mbh > 10^ Mq) QSOs at 
z~6, leading to a discussion that super massive BHs at z^b grow 
rapidly without commensurate growth of their host galaxies. On 
the contrary, this QSO is consistent with the local M-a relation as 

shown in Fig[5] 

Based on their galaxy merger simulations, iRobettson et al.l 

( 12006 ^ predicted a weak redshift-dependent shift in the M-cr re- 
lation due to an increasing velocity dispersion for a given galactic 
stellar mass. This QSO is even consistent with their predicted M- 
cr relation at z=6 shown in the purple dotted line in Fig[5] Note, 
however, that our measurement is made under the assumption that 
the gas is virialized. It is assumed that there is no effect on the 
measured line profile from non-kinematic broadening mechanisms 
(such as resonant atomic scattering or dust scattering). 

We do not claim an evolution or non-evolution of the 
MBn/Mbuige ratio based on one data point at z=6.4. However, 
due to the optical selection limit, previous QSO samples at z~6 
were limited to extremely massive BHs with Mbh > W'^Mq, 
sampling only a small range in BH mass. Combined with errors 
in measuring velocity dispersion, this may have been one of the 
reasons why an almost flat Mbh — cr relation was found at z~6 
jShields et al.l2006l : IWang et alllOlG ). and shown in FigO Our re- 
sult with a fainter QSO highlights the importance of sampling a 
wider range in Mbh to more accurately assess the Mbh — o" re- 
lation at z~6. Extended Lya emission around a QSO may offer an 
alternative way to investigate fainter QSOs, whose molecular gas 
lines cannot be easily observed with existing facilities. 



4.4 Note added at revision 

After this paper was submitted and a referee's report received, an 
observati onal spectroscopic s tudy of the identical QSO appeared 
on arXiv. Iwillott et alj I20T If) observed the same Ly-a halo with 
Keck/ESI for 8.5 hours in seeing of 0.93" They found hLya > 
8 X 10*^ erg/s, with an extension of over 15 kpc, and a FWHM 
of 640 km/s. However, in contrast to the work presented here, an 
independent PSF was not available to remove the underlying QSO 
continuum. We show here that the results are very similar if a PSF 
star or the QSO continuum are utilised to remove the QSO across 
the Ly-Qf line, thus independently strengthening the result in Willott 
et al. In consequence the numbers for the extended Ly-a flux and 
width reported in Willott are consistent with the findings in this 
paper. 

We thank the referee and M.Koss for insightful comments. 
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